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Abstract. In recent years, with the advent of the web of data, a growing num¬ 
ber of national mapping agencies tend to publish their geospatial data as Linked 
Data. However, differences between traditional GIS data models and Linked Data 
model can make the publication process more complicated. Besides, it may re¬ 
quire, to be done, the setting of several parameters and some expertise in the se¬ 
mantic web technologies. In addition, the use of standards like GeoSPARQL (or 
ad hoc predicates) is mandatory to perform spatial queries on published geospa¬ 
tial data. In this paper, we present GeomRDF, a tool that helps users to convert 
spatial data from traditional GIS formats to RDF model easily. It generates ge¬ 
ometries represented as GeoSPARQL WKT literal but also as structured geome¬ 
tries that can be exploited by using only the RDF query language, SPARQL. Ge¬ 
omRDF was implemented as a module in the RDF publication platform Datalift. 
A validation of GeomRDF has been realized against the French administrative 
units dataset (provided by IGN France). 


1 Introduction 

Over the last decade, efforts to share geospatial data on the Web have mainly focused 
on the development and the standardization by ISO TC 211 and the Open Geospatial 
Consortium (OGC) of domain-specific good practices and tools called Spatial Data In¬ 
frastructures (SDIs) [5]. These good practices have been particularly confirmed by the 
adoption by the European Parliament of the INSPIRE Directive establishing an infras¬ 
tructure for spatial information in the Community [22]. At the same time, generic good 
practices for data sharing on the Web called Linked Data have been developed and 
supported by the W3C [3]. Over the last years, many datasets have been published ac¬ 
cording to the Linked Data principles, so that these good practices have spread and are 
now being adopted for geospatial data. Eor example, Geo.LinkedData.es [10], is an ini¬ 
tiative to enrich the Web of data with Spanish geospatial data. Besides, LinkedGeoData 
[18], is an effort to add information collected by the OpenStreetMap project to the Web 
of data. Moreover, the Ordnance Survey Linked Data [21] is an initiative of the Great 
Britain’s national mapping agency to publish a number of its products as Linked Data. 

Although these initiatives are promising, the amount of geospatial data published as 
Linked Data remains limited compared to available geospatial data published through 
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other means like Spatial Data Infrastructures or Open Data platforms. This could be 
partly due to the fact that converting geospatial data from their traditional GIS formats 
to fully interconnected RDF data raises a domain-specific difficulty: converting into 
RDF the geometries used by geographic information standards for representing the lo¬ 
cation and the shape of real world geographic entities. This difficulty occurs at the first 
step of the publishing process. To be overcome, it requires not only geolocation vocab¬ 
ularies, but also specific tools that automatise the conversion of existing geospatial data 
into Linked Data [31], 

To deal with these issues, we present the GeomRDF tool which offers to geospatial 
data publishers the possibility to publish (in an easy-way) their data, and to users the 
possibility to query these geospatial data, using semantic web technologies. Most im¬ 
portantly, the GeomRDF tool bundles a fine-grained representation of geometry based 
on a vocabulary [11] that re-uses and extends the existing geographic standards and 
vocabularies (GeoSPARQL [13], NeoGeo [19]). This representation enables automated 
agents to reason over geometries. GeomRDF is implemented as a module^ in the Datalift 
Platform [7] in order to provide a complete path from geospatial data (ESRI Shape- 
file,DBMS and GML) to fully interlinked, identified, and qualified linked data. 

The paper is organized as follows. In the next section, we present some related 
works. In section 3 we present in detail the main components of the GeomRDF tools, 
together with the supported geospatial formats and external libraries. In section 4 we 
present a test of these tools on a dataset that represents French administrative units, 
provided by IGN France. Finally we conclude and give some perspectives in Section 5. 

2 Related Works 

Real world phenomenon are represented in geographic vector databases by geographic 
features, described by thematic properties and a geometry. Thus, converting geospatial 
data from their original standard to RDF implies representing not only their thematic 
properties - which can be handled like any other data, but also their associated geome¬ 
tries in RDF. Many knowledge extraction and data conversion tools have been proposed 
to generate RDF data from unstructured or structured data sources [26] [4] [6] [1] [29] 
[24] [17], but just a few vocabularies and their associated conversion tools have been 
implemented to address the issue of converting easily geospatial data into the RDF 
model. 

Geometry2RDF [12] is a tool, developed in the context of Geo.LinkedData.es [10], that 
enables geospatial data conversion into RDF. The conversion process takes as input 
different geospatial formats (ESRI Shapefile, GML and geospatial DBMS) and gen¬ 
erates RDE triples according to the NeoGeo vocabulary [19]. Consistently with Neo- 
Geo vocabulary, geometries associated to geographic features are therefore represented 

^ The current release of the Datalift Platform (available at: https : / /gf orge . inria . f r/ 
scm/?group_id=2935) includes the experimental version of GeomRDF (that does not 
contain the geometry structuring process). The final version of GeomRDF code source will be 
published with the next Datalift release. 
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in RDF as structured geometries. NeoGeo vocabulary defines geometries following 
OGC’s GML Simple Features Profile. This vocabulary directly reuses the class Point of 
the W3C Geo Vocabulary [30], and its related properties. As a consequence, geometries 
represented with NeoGeo vocabulary can be defined only with WGS84 coordinates. 

TripleGeo [23] is an extension of Geometry2RDF that takes into account, in the 
transformation process, the GeoSPARQL vocabulary. This vocabulary for represent¬ 
ing RDF geospatial data has been defined as part of the OGC standard GeoSPARQL, 
which also provides an extension to the SPARQL query language for processing RDF 
geospatial data. This vocabulary defines two RDFS datatypes , namely http: / /www. 
opengis.net/ont/geosparql#wktLiteral and http://www.opengis. 
net/ont/geosparql#gmlLiteral, for representing geometries respectively as 
WKT serialization as defined by the Simple Feature standard [16] and as GML serial¬ 
ization as defined by GML standard [20]. Both serializations imply expliciting for each 
geometry the coordinate reference system (CRS) used. For that purpose, OGC main¬ 
tains a set of CRS URIs for identifying the most common CRSs. In addition, TripleGeo 
provides an on-the-fly functionality for transforming geometries coordinates from one 
given coordinate reference system to another. 

shp2GeoSPARQL [27] is also an extension of Geometry2RDF which transforms ge¬ 
ometrical information from geospatial datasets to RDF. shp2GeoSPARQL parses ESRI 
Shapefile in order to retrieve the geometries of features, and generates a RDF represen¬ 
tation of geometries consistent with the GeoSPARQL geometry vocabulary. 

These three solutions have their advantages and drawbacks. On the one hand, the 
tools based on GeoSPARQL vocabulary have the advantage of allowing geometries de¬ 
fined in any CRS. However, there are very few triplestores that implement this standard 
yet [2], so that GeoSPARQL spatial operators and functions can not be used commonly 
with the available linked geospatial data. On the other hand, the tool based on NeoGeo 
vocabulary has the advantage of providing structured geometries that can be handled 
with regular SPARQL queries. But geometries coordinates are restricted to WGS84 
CRS. In order to overcome these limitations, we propose the GeomRDF tool. It is based 
on a vocabulary that reuses and extends GeoSPARQL simple feature vocabulary and 
NeoGeo so that geometries can be defined in any CRS and represented both as struc¬ 
tured geometries that can be handled with regular SPARQL queries and as wktLiteral 
datatypes that are compliant with GeoSPARQL standard. 


3 GeomRDF Components 

In this section we describe the different components of GeomRDF (cf. Fig. 1) and 
the process that transforms a geospatial dataset from traditional GIS formats to RDF. 
GeomRDF includes three components: the input parser, the features parser and the RDF 
builder. We outline in the following sections these components. 


3.1 The Input Parser: 

GeomRDF takes as input ESRI Shapefile [8], Geospatial DBMS or GML [20]. This 
component extracts from each file/database all the features and their descriptions (e.g.. 
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schema, types, CRS (Coordinate Reference System). We mean by ’’feature”, the ’’Sim¬ 
ple Feature” defined by the Open Geospatial Consortium (OGC) [25] and the Interna¬ 
tional Organization for Standardization (ISO) standard ISO 19125 to have both spatial 
and non-spatial attributes. Spatial attributes are geometry valued, and simple features 
are based on 2D geometry with linear interpolation between vertices (e.g., Multipoly- 
gones. Polygons, .etc). 


ESRI Shapefile: This format is a part of the GIS technology offered by ESRI. Spatial 
data format defined by ESRI stores non-topological geometry and attribute information 
for the spatial features in a dataset. 

Regarding this format, four hies are processed: the three mandatory dBASE hie (.dbf), 
index hie (.shx) and main hie (.shp), plus the metadata hie (.prj) that describes the CRS 
used by the dataset (In the case where this hie is missing, the default considered CRS is 
WGS84). The schema and the type of the thematic and geometric attribute are extracted 
from the main and dBase hie. 


Geospatial DBMS: It is a kind of database that is optimized to store and query geospa¬ 
tial data. GeomRDF can handle: IBM DB2, H2, MySQL, Oracle Spatial, PostGIS (a 
spatial extension to PostgreSQL), SpatiaLite (a spatial extension to SQLite), Microsoft 
SQL Server and Teradata. 
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GML: OGC standard encoding specification for geodata in XML that enables the stor¬ 
age, transport, processing, and transformation of geographic information. 

These three format of storage are well known in the geographic information science 
field and are used by the majority of mapping agencies. 

3.2 Feature Parser: 

This component iterates on the features extracted during the first step and classifies 
their properties depending on their type, either thematic or geometric. For one feature 
the parser tests all the properties and, for each property, stores information like its name, 
its value and its type. When the property is geometric, other information concerning this 
property (e.g., the number of geometries contained by a Multipolygon) will be stored. 
The Fig. 2 shows the two classes designed to store feature properties. 



Fig. 2: Feature Properties Class Diagram 


To illustrate how features are processed, we take as example the department of 
’’Paris” defined in GEOFLA dataset on French administrative units produced by IGN 
France and stored in ESRl Shapefile format. We detail thereafter the different steps for 
storing thematic and geometric properties of this feature. 

Thematic properties All the thematic properties that describe the department of’’Paris” 
are stored as a ’’EeatureProperty” structure (cf. fig 2). These properties include, for in¬ 
stance, the department code (CODEJDEPT) the department name (NOM_DEPT), the 
region code (CODE_REG) and the region name (NOM_REGION) (the other properties 
are detailed in the specifications of this dataset"^). 

http://professionnels.ign.fr/sites/default/files/DC_GEOFLA_l- 
1. pdf 












6 


Fay§al Hamdi, Nathalie Abadie, Benedicte Bucher, and Abdelfettah Feliachi 


The following table shows a sample of the stored thematic properties, their values and 
their types: 


Property 

CODE_DEPT 

NOM_DEPT 


CODE_REG 

NOM_REGION 

Value 

75 

PARIS 


11 

ILE DE FRANCE 

Type 

String 

String 


String 

String 


Table 1: The thematic properties of the department of ’’Paris” 


Geometric properties The geometry associated to the department of ’’Paris” is stored 
as a sequence of coordinates defining the points that are part of the boundary of the 
department. Because of the length of this geometry, we represent only some points: 

MULTIPOLYGON(( (2.41633 48.84923, 2.41597 48.84662, .... , 2.41633 

48.84923))) 

This geometry is stored as ’’GeometryProperty” structure (cf fig 2). In this exam¬ 
ple, the geometry of the department of’’Paris” is represented by a Multipolygon. Indeed, 
this type of geometry is a good example because it includes the majority of the other 
types of geometry handled by GeomRDF. In this case, the parser stores at first all the 
Polygons that compose the original Multipolygon. Then, it stores the exterior and the 
eventual interior LinearRings that compose each polygon, as well as all the points in¬ 
cluded in each LinearRing. Finally the coordinates of each point are stored. 

At this stage, GeomRDF provides a very useful functionality which consists in 
transforming these coordinates from their original CRS (Coordinate Reference Sys¬ 
tem) to another one. This transformation is performed on-the-fly; before storing the 
geometry, the parser check whether there is a need or not to convert the CRS. 

3.3 RDF Builder: 

The RDF Builder generates, from the parsed features, a collection of RDF triples ex¬ 
pressed as subjects, predicates, and objects. As in the Feature Parser component, the¬ 
matic and geometric properties are treated separately. We illustrate in the following, the 
process of generating RDF triples from the feature (the department of ’’Paris”), parsed 
in the above example. 

Thematic properties Consider the following base URI data: <http: //data. ign 
. fr/id/geof la/departement/>. The value of a property is captured as text with 
the appropriate RDF literal type. For instance the value of the NOM/REGION will be 
represented by the literal " ILE de France" ~ "xsd: string. 

By default, GeomRDF generates predicates by reusing properties names. They can be 
replaced afterwards by predicates from the dataset vocabulary. This task of matching 
and replacing default predicates with well defined ontology predicates is handled within 
another module of the Datalift platform. Listing 1.1 shows an example of RDF triples 
representing thematic properties. 
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@prefix rdf: <http://www.w3.Org/1999/02/22-rdf-syntax-ns#> 

@prefix geofla: <http://data.ign.fr/def/geofla#> 

<http://data.ign.fr/id/geofla/departement/75> a geofla:Departement . 

<http://data.ign.fr/id/geofla/departement/75> rdfs:label "PARIS"@fr . 

<http: / /data. ign. f r/id/geof la/departeinent/75> geof la: codeDpt "75" '' 'xad: string . 


Listing 1.1: A sample of RDF thematic triples of the department of ’’Paris” 

Geometric properties The generation of the geometric RDF triples is guided by a vo¬ 
cabulary [11] that extends GeoSPARQL standard and gives a very precise description 
of the different types of geometry. Based on this vocabulary, GeomRDF generates, for 
each geometric property stored by the Feature parser, the appropriate RDF triples that 
represent the geometry and its CRS. Furthermore, as the the vocabulary of geometries 
was built in OWL language, a semantic reasoner can check the consistency and infer 
new axioms against generated triples. Each instance of a geometric class is linked to the 
URL identifying the CRS used for its coordinates thanks to the geom: cr s predicate. 
To well detail the functionality of this component we take as example the process of 
generating RDF triples from the geometric property of the department of ’’Paris”. 


In order to structure a Multipolygon, we follow the axioms defined in the vocabu¬ 
lary of geometries. We outline thereafter, the axioms used to define a Multipolygon and 
some samples of the generated RDF triples. 


Axiom 1 (Multipolygon) An instance of the class geom:MultiPolygon is as¬ 
sociated to one or more other instances of the class geom: Polygon via the property 

geom: polygonMember. 

@prefix geom: <http://data.ign.fr/def/geometrie#> 

<http://data.ign.fr/id/geofla/departement/75> a geom:MultiPolygon ; 
geom:crs <http://data.ign.fr/id/ignf/crs/WGS84GDD> ; 
geom:polygonMember _:bn00001 ; 
geom:polygonMember _:bn00002 , 


Listing 1.2: A sample of RDF triples representing a Multipolygon 

Axiom 2 (Polygon) An instance of the class geom:Polygon is defined by its 
contour which is an instance of the class geom: LinearRing associated with the in¬ 
stance of geom: Polygon via the property geom: exterior. A polygon has a single 
contour. It may also contain holes (0 or more), defined as instances of 

geom: LinearRing associated to geom: Polygon via the property 

geom: interior. 

_:bn00001 a geom:Polygon ; 

geom:crs <http://data.ign.fr/id/ignf/crs/WGS84GDD> ; 
geom:exterior _:bn000021 . 

_:bn00002 a geom:Polygon ; 

geom:crs <http://data.ign.fr/id/ignf/crs/WGS84GDD> ; 
geom:exterior _:bn000031 . 


Listing 1.3: A sample of RDF triples representing Polygons 
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Axiom 3 (LinearRing) An instance of the class geom: LinearRing is associ¬ 
ated to an instance of geom: Point sList, which is a list of points, via the prop¬ 
erty geom: points. This instance of PointsList must be associated to its initial point 
(which is also its last point) via the property geom: firstAndLast. 

_:bn000021 a geom:LinearRing; 

geom:crs <http://data.ign.fr/ici/ignf/crs/WGS84GDD> ; 
geom:points _:bn000022 ; 
geom:firstAndLast _:bn000023 . 

_:bn000031 a geom:LinearRing; 

geom:crs <http://data.ign.fr/id/ignf/crs/WGS84GDD> ; 
geom:points _:bn000032 ; 
geom:firstAndLast _:bn000033 . 


Listing 1.4; A sample of RDF triples representing LinearRings 

Axiom 4 (PointsList) A geom:PointsList is a subclass ofrdfiList. An instance 
of geom: Point sList is composed of only instances of type geom: Point. 

_:bn000022 a geom:PointsList; geom:firstAndLast _:bn000023 ; 
rdf:rest _:bn000213 . 

_:bn000213 a geom:PointsList; rdf:first _:bn000024 ; 
rdf:rest _:bn0000214 . 

_:bn000219 a geom:PointsList; rdf:first _:bn000210 ; 
rdf:rest _:bn000022 . 


Listing 1.5: A sample of RDF triples representing a PointsList 

Axiom 5 (Point) An instance of geom: Point is associated to its coordinates via 
the properties geom: coordX and geom: coordY. 

_:bn000023 a geom:Point ; 

geom:crs <http://data.ign.fr/id/ignf/crs/WGS84GDD> ; 
geom:coordX "2.41633”"'xsd:double ; 
geom: coordY ”48.84923"'' "'xsd: double . 


Listing 1.6: A sample of RDF triples representing a Point 

In addition to this structured representation of geometry coordinates, the geome¬ 
try will be associate to a WKT literal according to the GeoSPARQL standard, as the 
following: 

@prefix geosparql: <http://www.opengis.net/ont/geosparql#> 

<http://data.ign.fr/id/geofla/departement/75> geosparql:asWKT "<http://data.ign.fr 
/id/ignf/crs/WGS84GDD> MULTIPOLYGON({(2.41633 48.84923, 2.41597 48.84662, 

.... , 2.41633 48.84923)))”“ “<http://www.opengis.net/ont/geosparql#wktLiteral 


Listing 1.7: A sample of an RDF triple representing WKT literal 


3.4 External Libraries 

The actual version of GeomRDF (a platform-independent version implemented in Java) 
depends on only two essential libraries; GeoTools [15] and OpenRDF Sesame [14]. 
Their descriptions are listed below: 
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GeoTools is an open source Java library that provides methods for manipulating geospa¬ 
tial data. It includes many features such as a support for a wide number of formats 
(including ESRI Shapefile, XML (GML, KML), and geospatial databases), and a co¬ 
ordinate transformation service, providing coordinate reference systems descriptions 
and operations on coordinates such as map projections and datum transformations. The 
data structures used in GeoTools are based on Open Geospatial Consortium (OGC) 
standards. 


OpenRDF Sesame is an open source Java framework for processing RDF data. Its fea¬ 
tures includes parsers, storage solutions (Triplestores), reasoning and querying, using 
the SPARQL query language. Sesame supports all main RDF file formats, including 
RDF/XML, Turtle, N-Triples, TriG and TriX. 

In an experimental version of GeomRDF we have tested some command-line avail¬ 
able in GDAL/OGR [9] library. We noticed that these command-line makes our tool 
dependent on native library, and hence less portable. Therefore, we have decided to rule 
out the use of this library. 


4 Transforming the French Administrative Units Dataset from 
GIS Format to RDF Model 


This section presents a test of the conversion and the publication of the French Ad¬ 
ministrative Units dataset, using the GeomRDF tool and the Datalift platform. This 
dataset, namely GFOFLA, is provided by the IGN France as a set of FSRI Shapefiles. 
It contains the description of all the administrative subdivisions (department, commune, 
arrondissement and canton) of France, overseas departments and Mayotte. 

To convert and publish GFOFLA we follow the following steps: (note that all these 
steps are realized within the Datalift platform). 


Conversion of geospatial data to RDF In this step we used the GeomRDF module to 
convert the FSRI Shapefiles to RDF. A transformation from the original CRS of the ge¬ 
ometries (Lambert-93) to WGS84 is carried out. This choice is motivated by the fact that 
most of the available geospatial linked datasets use this CRS. For identifying resources, 
a base URI is constructed (according to W3C best practice) as the following: <http 
://data.ign.fr/id/geofla/[types_of_administrative_division]/[id]>. For example the com¬ 
mune of ’’Saint-Mande” will be identified by: <http://data.ign.fr/id/geofla/commune 
/94067>. 

In order to illustrate the fact that geometries can be described with different CRSs, 
the points representing respectively the location of the head office and the location of 
the centroid of administrative units are represented by Lambert-93 coordinates. See for 
example the location of the head office of the ’’Seine et Marne” department: http: 
//data.ign.fr/id/geofla/departement/Point_94 02 8. 
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Mapping RDF data to vocabularies An ontology describing the GEOFLA dataset 
was developed and published by the IGN. Two modules of the Datalift platform can 
be used to map the RDF data, generated in the previous step, to the ontology classes 
and properties: ”RDF conversion based on ontologies” or ’’RDF2RDF transformation”. 
We have used the latter which performs SPARQF CONSTRUCT queries in order to 
generate a new data named graph according to a given pattern. Fisting 1.1 (section 3) 
shows an example of RDF triples transforming according to the GEOFFA ontology. 

Publication of tbe datasets In this step, the RDF data are published online as Finked 
Open Data and accessible through the following SPARQF endpoint address: http: 
//data.ign.fr/id/sparql. 


Interlinking with other datasets As shown in the listing 1.9, we have used the ’’RDF2RDF 
transformation” module and administratives units identifiers in order to add owksameAs 
links with the INSEE (French National Institute for Statistics and Economic Research) 
web dataset. A data interlinking module based on the Silk Framework [28] is also avail¬ 
able within the Datalift platform when more complicated interlinking processes are 
needed. 

@prefix owl: <http://www.w3.Org/2002/07/owl#> 

<http://data.ign.fr/id/geofla/arrondissement/751> owl:sameAs <http://id.insee.fr/ 
geo/arrondisseinent/751> 

<http://data.ign.fr/id/geofla/departeinent/75> owl:sameAs <http://id.insee.fr/geo/ 
departement/75> 


Fisting 1.8: sameAs links between the administrative unit of Paris in IGN and INSEE 

Spatial queries over structured geometries As explained in section 2, we have de¬ 
cided to generate structured geometries in order to enable users perform basic spatial 
queries in SPARQF, even without storing our geospatial data in a triplestore that imple¬ 
ments GeoSPARQF standard. The following query illustrates that scenario. It returns 
all the instances of rdPType geofla:Departement whose geometries intersect a bounding 
box whose coordinates are given as parameters of the query. 

PREFIX rdf: <http://www.w3.org/1999/02/22-rdf-syntax-ns#> 

PREFIX rdfs: <http://www.w3.org/2000/01/rdf-schema#> 

PREFIX geofla: <http://data.ign.fr/def/geofla#> 

PREFIX geom: <http://data.ign.fr/def/geometrie#> 

SELECT DISTINCT ?name WHERE { 

?dep rdf:type geofla:Departement . 

?dep rdfs:label ?name . 

?dep geom:geometry/geom:polygonMember/geom:exterior/geom:points ?pl . 

?pl rdf:type geom:PointsList . 

?pl (rdf:rest*/rdf:first)|geom:firstAndLast ?pm . 

?pm rdfrtype geom:Point . 

?pm geom:coordX ?x . 

?pm geom:coordY ?y . 

FILTER ((?x > 1 && ?x < 3) &&(?y > 42 && ?y < 44)) . 

} 


Fisting 1.9: Example of a spatial query over structured geometries in SPARQF 
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This query can be applied directly on GEOFLA dataset thanks to this SPARQL 
Endpoint http : //data . ign . fr/id/sparql. Results are the departments called 
’’Ariege”, ’’Tarn”, ”Tarn-et-Garonne”, ”Aude”, ’’Aveyron”, ’’Pyrenees-Orientales”, ’’Haute- 
Garonne”, ”Gers” and ’’Herault”. 


5 Conclusion and Future Work 

In this paper we presented the GeomRDF tools that enables the transformation of 
geospatial data into RDF. We presented the input format that GeomRDF supports, its 
mains components, and its dependencies on external libraries. We detailed how the ge¬ 
ometric informations are structured according to an ontology that re-uses and extends 
the existing geographic vocabularies, and showed an example of this structuration in 
the case of a Multipolygon geometry. 

The GeomRDF tool has been implemented as a module in the Datalift platform. We 
illustrated an application of this tool by detailing the process of publishing the French 
Administrative Units dataset provided by IGN France. 

So far GeomRDF takes only geospatial formats as input data. However, it would 
be useful to make it also accept text file formats, such as CSV, containing direct loca¬ 
tion information described by coordinates, and convert these textual coordinates into 
instances of the class geom:Point. Then, such location information could be handled by 
query, interlinking or visualization tools in the same way as any other piece of geospa¬ 
tial information. 
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